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Abstract Infections by type I feline coronaviruses 
(FCoVs) have been shown to be significantly correlated 
with fatal feline infectious peritonitis (FIP). Despite nearly 
six decades having passed since its first emergence, dif- 
ferent studies have shown that type II FCoV represents 
only a small portion of the total FCoV seropositivity in 
cats; hence, there is very limited knowledge of the evolu- 
tion of type I] FCoV. To elucidate the correlation between 
viral emergence and FIP, a local isolate (NTU156) that was 
derived from a FIP cat was analyzed along with other 
worldwide strains. Containing an in-frame deletion of 442 
nucleotides in open reading frame 3c, the complete genome 
size of NTU156 (28,897 nucleotides) appears to be the 
smallest among the known type II feline coronaviruses. 
Bootscan analysis revealed that NTU156 evolved from 
two crossover events between type I FCoV and canine 
coronavirus, with recombination sites located in the 
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RNA-dependent RNA polymerase and M genes. With an 
exchange of nearly one-third of the genome with other 
members of alphacoronaviruses, the new emerging virus 
could gain new antigenicity, posing a threat to cats that 
either have been infected with a type I virus before or never 
have been infected with FCoV. 


Keywords Type II feline coronavirus - Feline infectious 
peritonitis - Genome analysis - Virus recombination - Virus 
evolution 


Introduction 


Feline coronaviruses (FCoVs) are large, enveloped, posi- 
tive-strand RNA viruses with a genome of approximately 
29,200 nucleotides [1-3]. The FCoVs belong to the genus 
Alphacoronavirus, family Coronaviridae, order Nidovi- 
rales. Other members of this subgroup include canine 
coronavirus (CCoV), transmissible gastroenteritis virus 
(TGEV), raccoon dog CoV (RDCoV/GZ43/03), and Chi- 
nese ferret badger CoV (CFBCoV/DM95/03) [4]. 

FCoVs are associated with diseases that range from sub- 
clinical and/or mild enteric infections to fatal infectious 
peritonitis [5]. Despite the high prevalence of FCoVs in feline 
populations around the world, only 5-12 % of seropositive 
cats develop feline infectious peritonitis (FIP). FIP is a 
chronic, progressive, immune-mediated disease in domestic 
and nondomestic fields. The typical histopathological finding 
of this disease is systemic perivascular necrotizing pyogra- 
nulomatous inflammation [6]. Two serotypes that differ in 
their growth characteristics in tissue culture and in their 
genetic relationship to CCoV and TGEV have been identified 
[7, 8]. Type If FCoV is significantly correlated with FIP when 
compared to type I viruses [9, 10]. However, unlike the 
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ubiquity of type I FCoV, infection by type II virus encom- 
passes only a small percentage of the total number of FCoV- 
seropositive cats in different studies [9-16]. 

Type II FCoV is estimated to have diverged from 
alphacoronavirues in 1953 [4]. Based on partial genomic 
sequence analysis, type II FCoVs were suggested to result 
from a double recombination between type I FCoV and 
CCoV [17]. Like most RNA viruses, CoVs mutate at a high 
frequency due to the high error rate of RNA polymeriza- 
tion. In addition, a unique feature of CoV genetics is the 
high frequency of RNA recombination in the natural evo- 
lution of this virus [18]. Recombination among CoVs is an 
attribute of the genus and is thought to contribute to the 
emergence of new pathotypes, such as severe acute respi- 
ratory syndrome CoV [19, 20], human CoV NL63 (HCoV 
NL63) [21], HCoV HKUI1 [22], and avian infectious 
bronchitis virus IBV) [23-26]. To gain better evolutionary 
insight into type II FCoVs, we analyzed the complete 
genome of a novel type II FCoV isolate. Taking our data 
together with data from other strains, we discuss the evo- 
lution of type II FCoV. 


Materials and methods 
Virus and isolation of viral RNA 


FCoV NTUI56 was isolated in 2007 from a kitten with 
naturally occurring FIP by the co-cultivation of pleural 
effusion with feline fewf-4 cells [27]. After three rounds of 
purification by limited dilution, the virus was propagated 
and titrated. All of the viruses used in this study for the 
sequencing of complete genome came from a stock virus 
passaged six times. NTUI56 is relatively fast-growing, 
induces a coronavirus-typical syncytial cytopathic effect 
and is a type II FCoV [10]. 


Reverse transcription and PCR amplification of viral 
sequences 


Eleven microliters of isolated RNA was added to the pre- 
mix, consisting of 4 ul of Sx RT buffer, 2.5 mmol dNTPs 
(GeneTeks Bioscience, Inc., Taipei), 50 pmol random 
primer, 0.2 mol dithiothreitol, and 1 ul of 200 U Moloney 
murine leukemia virus reverse transcriptase (Invitrogen, 
CA, USA) in a 0.6 ml reaction tube. This reaction mixture 
was then briefly centrifuged and incubated at 37 °C for 
60 min, then at 72 °C for 15 min, and finally at 94 °C for 
5 min. A total of 120 primers for PCR were chosen from a 
relatively conserved region of the FCoV genome. Follow- 
ing reverse transcription, | pl of the RT reaction mixture 
was added to 49 ul of the PCR mixture, which consisted of 
5 ul of 10x Tag buffer, each primer (10 pmol), dNTP 
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(2.5 mmol), 2 U of Tag DNA polymerase (GeneTeks, 
BioScience, Inc., Taipei), and 39 wl of 0.1 % DEPC water. 
An ABI-2720 thermal cycler (Applied Biosystems, USA) 
consisted of 3 min of preheating at 94 °C, followed by 35 
cycles of denaturation at 94 °C for 30 s, annealing at 55 °C 
for 30s, and extension at 72 °C for 1 min with a final 
extension at 72 °C for 7 min. The viral RNA termini were 
amplified using 3’- and 5’-rapid amplification of cDNA 
ends (Invitrogen, USA). 


Analysis of PCR-amplified products and sequencing 


A total of 10 ul of PCR products from each PCR mixture 
was analyzed using a 1 % agarose gel (Viogene, Taipei) for 
electrophoresis. Amplification products were visualized 
using UV illumination after ethidium bromide staining. 
The nucleotide sequences of the targeted DNA fragments 
were purified (Geneaid Biotech, Taipei) and sequenced in 
both directions using an auto sequencer (ABI 3730XL, 
USA). Full-length genome sequencing of NTU156 was 
performed by single-round PCR with a set of overlapping 
PCR products (average size 750 bp) that encompassed the 
entire genome. 


Sequence analysis 


The complete sequences of NTU156 were then compared 
with other alphacoronaviruses and the results are sum- 
marized in Table 1. Multiple alignments of nucleic acid 
sequences were performed by the Clustal W method using 
the MegAlign program (DNASTAR Inc., WI, USA). 
Phylogenetic analyses were conducted using MEGA, ver- 
sion 4.0. Similarity graphs were prepared with SimPlot 2.5 
software [28]. Potential recombination sites were identified 
using the Recombination Detection Program (RDP) [29]. 


Results 
Genomic sequence of FCoV NTU156 


The full genomic RNA sequence of FCoV NTU156 com- 
prises 28,897 nucleotides (nts), excluding the 3’ poly- 
adenylated nts. Sequence analysis revealed that NTUI56 
contains conserved open reading frames with an overall 
genome organization similar to known FCoVs (Table 2). 
The overall nucleotide composition is as follows: A, 
29.3 %; C, 17.3 %; G, 21.0 %; and T, 32.4 %. The G+C 
content 1s 38.3 %. NTUI56 possesses the putative tran- 
scription regulatory sequence (TRS) motif, 5’-CUAAAC-3’, 
at the 3’ end of the leader sequence and preceding each ORF 
(Table 2). 
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Table 1 Isolation dates and accession numbers of NTU156 and other alphacoronavirses 


Isolation 
Virus Host Year Country 
FCoV NTU156°* Cat 2007 Taiwan 
FCoV NTU2 Cat 2003 Taiwan 
FCoV 79-1146" Cat 1979 USA 
FCoV Black* Cat 1980 USA 
FCoV UU2* Cat 1993 USA 
FCoV DF-2* Cat 1981 USA 
FCoV ClJe* Cat 2007 UK 
FCoV 79-1683" Cat 1979 USA 
FCoV KU-2 Cat 1991 Japan 
FCoV Yayoi Cat 1981 Japan 
FCoV C3663 Cat 1994 Japan 
CCoV NTU336° Dog 2008 Taiwan 
CCoV 1-71" Dog 1971 Germany 
CCoV TN-449* Dog 1980s USA 
CCoV 174/06 Dog 2006 Italy 
TGEV Purdue“ Pig 1946 USA 
TGEV M6* Pig 1965 USA 
TGEV TS* Pig 2004 China 
HCoV 229E* Human 1963 


“ Full genome available 
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GenBank accession no. 
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Table 2 Coding potentials and putative transcription regulatory sequences (TRS) of FCoV NTU156 


ORFs Start—end (nucleotide No. of nucleotides 
position) 
ORF la 313-2,417 12,105 4,035 
ORF lab 313-20,417 20,105 6,701 
(shift at 12,381) 
S 20,414—24,775 4,362 1,453 
ORF 3a 24,840—25,076 237 78 
ORF 3b 25,02 1—25,236 216 71 
ORF 3c 25233-25,433 201 66 
E 25,512—25,760 249 82 
M 25,77 1-26,556 786 261 
N 26,569-27,687 1,119 372 
ORF 7a 27,692-27,997 306 101 
ORF 7b 28,002—28,622 621 206 


5’, 3’ Untranslated regions (UTRs) and open reading 
frames (ORFs) 


The 5’-UTR of FCoV NTU156 comprises 312 nts—one or 
two nts more than other known FCoVs (79-1146: 311 nts, 
black: 311 nts, ClJe: 310 nts)—whereas the 3’-UTR 
of our virus comprises 275 nts, which is identical to other 
FCoVs. The 5’ two-thirds of the genome contain the la 
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No. of amino acids 


Putative TRS 


Nucleotide position TRS sequences 


94 GAACUAAACGAA210AUG 
20,382 UUACUAAACUUU23AUG 
24,832 GAACUAAACUUAUG 
25,469 GUUCUAAACGAA37AUG 
25,762 GAACUAAACAAAAUG 
26,557 UAACUAAACUUCUAAAUG 
27,684 GAACUAAACGCAUG 


(nt 313-12,417) and lab (nt 313—20,417) genes that encode 
the nonstructural polyproteins. A typical coronavirus “slip 
site,” 5’-UUUAAAC-3’ (nt 12,38 1—12,387) is located within 
this gene. These genes are followed by genes encoding the 
four structural proteins: spike (nt 20,414—24,775), envelope 
(nt 25,512—25,760), membrane (nt 25,771—26,556), and 
nucleocapsid (26,569—27,687) (Table 2). The two strings of 
accessory genes identified in all of the known FCoVs, 
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1.e., ORF 3ab and ORF 7ab, were found in NTU156 as well 
(Table 2). However, an in-frame deletion of 442 nucleotides 
in ORF 3c was identified, which resulted in a relatively short 
gene comprising only 201 nts. 


Genetic comparison and phylogenetic analysis 
with other alphacoronaviruses 


The overall sequence comparison revealed that NTUI56 
was more closely related to known subgroup la CoV but 
not Ib within alphacoronaviruses (Fig. 1). Nucleotide 
sequences similarity graphs of NTU156 with known type I 
FCoVs, CCoVs, and TGEVs were created by the Simplot 
software. The results showed that NTUI56 was more clo- 
sely related to type I FCoVs from the 5’ end of the genome 
to position 15,000 and from position 27,500 to 3’-UTR 
(Fig. 2). Genes located at the 5’ end (nsp1-11) and 3’ end 
(the N gene through ORF7) of NTU156 show consistently 
high similarity to type I FCoVs, whereas from nsp13 
through the E sequence, the similarity to canine and por- 
cine CoVs varies dramatically (Fig. 3). These data indicate 
that NTU156 might have arisen from recombination events 
between different strains of CoVs from species other than 
cats. Two possible recombination sites, at approximate 
positions 14,300 and 27,300, corresponding to the RNA- 
dependent RNA polymerase (RdRp) and the M gene, 
respectively (Fig. 4), were further analyzed. 

Phylogenetic trees using the nucleotide sequence of 
genes for putative proteins and polypeptides of alphacor- 
onaviruses were further constructed. At nsp 1 through nsp 
11 (Supplementary Fig. la) and from N gene (Supple- 
mentary Fig. 1b) to the ORF 7ab gene (Supplementary 
Fig. lc), NTU156 was more closely related to type I 
FCoVs. At the nsp 12 (RdRp) and the M gene, NTU156 
was not clustered with any known alphacoronavirues 
(Fig. 3a, f). From the nsp13 through the E gene, NTU156 


Fig. 1 Phylogenetic 
relationships constructed using 
the complete genome sequences 
of FCoV NTU156 and other 
alphacoronaviruses. Analysis 
was performed using MEGA 4 
software and neighbor-joining 
methods based on 1,000 
replicates. Bootstrap support 
values greater than 90 are 
shown 
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Fig. 2 Nucleotide sequence similarity of the complete genome of 
FCoV NTU156 with other alphacoronaviruses. A similarity plot was 
constructed to identify the sequence homology between type I FCoVs 
black, ClJe, and UU2 (gray); CCoV NTU336 (red); and TGEV 
purdue, M6, and TS (blue). Red arrows represent putative recombi- 
nation regions. A similarity of 1.0 indicates regions that share 100 % 
nucleotide identity. The similarity calculation was performed using 
the following parameters: a window size of 1,000 bp and a step size 
of 200 bp for the full-length sequence 


was clustered with CCoV (Fig. 3b-e). Taken together, 
these data indicated that NTU156 might have evolved from 
two recombination events with CCoV, with the sites of 
recombination located in the RdRp and M genes. 


Discussion 


A unique feature of CoV genetics is the high frequency of 
RNA recombination both in vivo and in vitro [18]. Here, an 
interspecies recombination between feline and canine CoV 
was identified in a viral strain NTU156, which was isolated 
from the pleural effusion of a FIP cat. This is the first time 
that evidence for natural recombination has been docu- 
mented through the complete genome sequence analysis of 
type II FCoV. In 1998, Herrewegh et al., based on partial 
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Fig. 3 Phylogenetic relationships constructed using nspl2 (a), 
nsp13-16 (b), S gene (c), ORF3 gene (d), E gene (e), and M gene 
(f) sequences of FCoV NTU156 and other alphacoronaviruses. 


sequence analysis, first determined that type II FCoVs 
79-1146 and 79-1683 originated from a homologous RNA 
recombination event between type I FCoV and CCoV [17]. 
The complete sequence of strain 79-1146 was later pub- 
lished in 2005 [1]. When comparing strain 79-1146, 
79-1683, DF-2, and our NTU156 strain, the only four type 
II FCoVs that have had their full-length genomes 
sequenced to date, a common phenomenon was found; both 
viruses arose from a double recombination event between 
type I FCoV and coronaviruses from other species. 
NTU156 appears to have evolved from a recombination 
between type I feline and canine CoV; however, when we 
aligned the genes in which recombination took place with 
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Analysis was performed using MEGA 4 software and neighbor- 
joining methods based on 1,000 replicates. Bootstrap support values 
greater than 90 are shown 


other type If FCoVs, genome crossovers with other 
alphacoronavirus were noted. When the sequences around 
the putative recombination sites were examined, 1.e., one 
located in the 5’ region (strain 79-1146) and two in the 3’ 
region (strain 79-1146, and 79-1683), porcine coronavirus 
(TGEV) was also found to contribute to the evolution of 
type II FCoV (in addition to CCoV) (Figs. 3b, 5). This 
finding is not surprising because the receptor for type I 
FCoV, feline aminopeptidase N, has been found to serve as 
a receptor for several CoVs, including canine, porcine, and 
human CoVs [30]. Therefore, the interspecies recombina- 
tion of type I FCoV with any of the above viruses might 
occur in nature. The recombination of type I FCoV in cats 
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Fig. 4 Bootscan analysis of the complete genome of FCoV NTU156 
with other alphacoronaviruses (type I FCoV UU2: pink; CCoV 
NTU336: yellow; TGEV purdue: green). Analysis was performed 
using RDP software and neighbor-joining methods based on 1,000 


living in the same household with dogs or living close to 
pig farms could give rise to type II FCoV. 

Based on the analysis of the four full genomes of type II 
FCoV available at present (NTU156, 79-1146, 79-1683, 
and DF-2), type Il] FCoVs appear to retain type I FCoV 
sequences in their 5’ and 3’ ends. We asked whether the 
genes located in these regions are indispensable for FCoV 
replication in cats. To answer that question, the amino acid 
sequence of genes retained at both ends of FCoVs, CCoVs, 
and TGEVs, 1.e., nsp 1 through nsp 11 and the N and ORF7 
genes, were further aligned and compared (Table 3). In 
contrast to the greater than 90 % of amino acid homology 
between different strains of FCoVs, the nsp 2, 3, and 6, as 
well as the N gene and the ORF7a gene of FCoVs, when 
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Fig. 5 Nucleotide sequence similarity of the complete genome of 
FCoV 79-1146 with other alphacoronaviruses. A similarity plot was 
constructed to identify the sequence homology between type I FCoVs 
black, ClJe, and UU2 (gray); CCoV NTU336 (red); and TGEV 
purdue, M6, and TS (blue). Red arrows represent putative recombi- 
nation regions. A similarity of 1.0 indicates regions that share 100 % 
nucleotide identity. The similarity calculation was performed using 
the following parameters: a window size of 1,000 bp and a step size 
of 200 bp for full-length sequences 


replicates. Putative recombination regions are indicated with arrow. 
The boostrap support calculation was undertaken using a window size 
of 1,000 bp and a step size of 200 bp. Dashed line indicates the 95 % 
bootstrap support threshold values 


compared with CCoVs or TGEVs, exhibited similarity of 
less than 80 %. This finding indicates that these gene 
products might possess irreplaceable functions for FCoV 
replication. This might explain why type II FCoVs found in 
nature harbor genomes that evolved from a double 
recombination. 

Although the prevalence of type II FCoV is consistently 
lower (2-11 %) than type I virus around the world 
(88-98 %) [9-16], our previous study indicates that 
infection of type II FCoV correlates significantly with FIP 


Table 3 Percent similarity with respect to the amino acid sequences 
of FCoV NTU156 and other alphacoronaviruses 


Gene FCoV CCoV TGEV 
NTU2 Black ClJe NTU336 174/06 Purdue TS 
ORF1 
nspl — 93.6 92.7 91.8 — 92.7 91.8 
nsp2. 9 — 93.6 91.7 80.2 — 79.1 78.7 
nsp3. — 91.5 84.9 79.2 — 78.8 78.3 
nsp¢  — 96.3 93.5 88.8 — 88.2 87.1 
nsp5  — 97.4 96.7 93.4 93.4 93.0 92.7 
nsp6 — 93.9 91.5 77.9 76.5 77.9 77.2 
nsp7 — 97.6 97.6 97.6 96.4 96.4 96.4 
nsp8 _— 97.9 96.4 92.8 90.8 92.8 92.3 
nsp9  — 95.5 96.4 86.5 91.0 89.2 89.2 
nsplO — 98.5 97.8 95.6 95.6 95.6 94.8 
nspll — 94.7 100 84.2 84.2 84.2 84.2 
N 93.3 92.0 91.7 77.5 77.2 76.4 76.9 
ORF7 
7a 99.0 97.1 96.1 81.4 81.4 78.5 74.7 
7b 88.9 86.9 86.7 59.9 61.8 — — 


* Not available 
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when compared to type I [10]. As shown in the present 
study, type II FCoV arises by exchanging a large genome 
fragment (approximately 12 kb) of type I FCoV with other 
members of alphacoronaviruses. The genes exchanged 
through this double recombination include nsp 13-16, 
structure protein S (spike), and accessory protein 3abc. The 
nsp 13-16 proteins are replication proteins with functions 
such as helicase activity (nsp 13), nucleoside triphospha- 
tase activity (nsp 13), RNA 5’-triphosphatase activity (nsp 
13), 3’-5’ exoribonuclease activity (nsp 14), RNA cap 
formation (nsp 14 and nsp 16), and endonuclease activity 
(nsp 15) [31]. It has ben reported that the function of the 3c 
protein might be crucial for viral replication in the gut but 
is dispensable for systemic FCoV replication [32]. How- 
ever, S proteins play a crucial role in receptor binding and 
eliciting protective immunity [33]. Through the replace- 
ment of nearly one-third of the genome, the new virus 
might gain new antigenicity, posing a threat to cats that 
either have been infected with a type I virus before or never 
have been infected with FCoV. 
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